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//-Stability of Nonlinear Positive Systems With 
Unbounded Time-Varying Delays 

Tianping Chen, Senior Member, IEEE, and Xiwei Liu, Member, IEEE 


Abstract —Stability of the zero solution plays an important role 
in the investigation of positive systems. In this note, we revisit 
the /x-stability of positive nonlinear systems with nnbounded 
time-varying delays. The system is modelled by continnous-time 
differential equations. Under some assumptions on the nonlinear 
fnnctions like homogeneous, cooperative, nondecreasing, we pro¬ 
pose a novel transform, by which the nonlinear system reduces 
to a new system. Thns, we analyze its dynamics, which can 
simplify the nonlinear homogenous functions with respect to 
(w.r.t.) arbitrary dilation map to those w.r.t. the standard dilation 
map. We finally get some criteria for the global /x-stability. 
A nnmerical example is given to demonstrate the validity of 
obtained results. 

Index Terms —Nonlinear, positive systems, time delay, un- 
bonnded, /x-stablllty. 

I. Introduction 

In many cases of the real world, states of systems should 
be always non-negative, for example, population levels in 
biology m, industrial processes involving chemical reactors 
El, transport and accumulation phenomena of substances in 
compartmental systems El-El. All these examples can be 
described mathematically by the model of positive systems 

E-El. 

In the literature of positive systems, the zero solution takes 
an important position, for example, the zero solution means 
the death of all species in biosphere. Moreover, although many 
important and interesting properties have been reported and 
analyzed, the most fundamental one is the stability property, 
which has been comprehensively studied by researchers from 
different fields. 

For a linear system to be positive x{t) = Ax{t), which is 
called positive linear time-invariant (LTI) system, a necessary 
and sufficient condition is that A should be a Metzler matrix 
(i.e., off-diagonal entries are nonnegative). For positive linear 
time-delay systems, El-El show that the presence of constant 
delays does not affect the stability performance of the system; 
Eol-IIIl (and ini-ini) also report this fact for positive 
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linear systems (and positive systems defined by functional 
and integro differential equations) of time-varying delays; 
m present some criteria for exponential stability of positive 
LTI differential systems with distributed delay; M gives an 
extension of the classical Perron-Frobenius theorem to positive 
quasi-polynomial matrices, then some necessary and sufficient 
conditions for the exponential stability of positive linear time- 
delay differential systems are obtained; in addresses the 
asymptotic stability of discrete-time positive systems with 
bounded time-varying delays and proves that the stability is 
also determined by the delay-free systems; establishes 
an equivalent relationship between asymptotical stability and 
exponential stability for discrete-time positive system for all 
bounded time-varying delays; mi studies the asymptotic sta¬ 
bility and decay rates with unbounded delays; and for positive 
linear switched system (PLS), mi-EQi analyze and prove that 
for 2-dimensional PLS, the necessary and sufficient condition 
for stability under arbitrary switching is that every matrix 
in the convex hull of the matrices defining the subsystems 
is Hurwitz, but is not true for 3-dimensional PLS; and ED 
also addresses the stability problem of both discrete-time and 
continuous-time PLSs with arbitrary (even unbounded) time 
delays. 

However, all these above papers study only the linear 
systems, which is just a small fraction of the whole positive 
systems. Recently, some great progress on the homogeneous 
systems fTlX has been made, which is a particular class of 
nonlinear systems. To overcome this disadvantage, there have 
been some excellent works on nonlinear positive systems in 
the literature. For example, ll23l is a pioneering work on 
generalizing the corresponding analysis from positive linear 
systems to homogeneous cooperative and irreducible systems; 
El shows that a constant-delayed homogeneous cooperative 
system is globally asymptotically stable (GAS) for all non¬ 
negative delays if and only if the system is GAS for zero delay; 
ESl shows that GAS and cooperative systems, homogeneous 
of any order with respect to arbitrary dilation maps are D- 
stable with constant time delay; flEj investigates the exponen¬ 
tial stability of homogeneous positive systems of degree one 
with bounded time-varying delays; and lIZTl generalizes the 
homogeneous positive systems to any degree, and the bounded 
time-varying delay to be unbounded, the asymptotic and p- 
stability are discussed for both continuous-time and discrete¬ 
time systems. 

In all, there are few papers considering the stability with 
unbounded time delay, see ED and IITTI . However, in our 
previous papers ll28l - E^ . we are the first to give explicitly 
the definition of /i-stability, which includes the exponential 
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Stability, the power-rate stability, the log stability, the log-log 
stability, etc. The key concept of defining a new Lyapunov 
function can even be retrieved in the first author’s paper 
lEO). Moreover, the proposed theory is successfully applied 
on the stability of equilibrium for neural networks En- 
(I33l, synchronization and consensus llTSl for complex 
networks. Therefore, in this note, we will revisit the /r- 
stability of nonlinear positive systems with unbounded time- 
varying delays using our proposed method, by defining a 
simple transform of the variables, we can directly get the 
corresponding criteria, which are easier to read and understand 
than previous works ll29l . 

The rest of this note is organized as follows. In section |II] 
some necessary definitions, lemmas and notations are given. 
In section [Till a simple transform is introduced and three 
properties are obtained. In section |IV] by transforming the 
original system to a simple form, whose nonlinear functions 
are homogeneous w.r.t. the standard dilation map, we analyze 
the ^-stability. Moreover, a numerical example to show the 
process of transform and analysis by using the obtained 
theoretical results is carefully demonstrated in Section [V] 
Finally, the note is concluded in section [Vl] 

II. Preliminaries 

Let R(R+),R" denote the field of real (positive) numbers 
and the nonnegative orthant of all n-dimensional real space 
R" respectively. If all elements of matrix A are non-negative 
(non-positive), then denote matrix A y 0(A 0). A real n x n 
matrix A = (a^), denoted as A G is Metzler if and 

only if its off-diagonal entries ^ j belong to M+. For 

two vectors x and y in R", if i = 1, • • • , n, where 

Xi iyi) is used to denote the i-th component of x (y), then we 
denote x > y. 

Definition 1: A dynamical system with state space R" is 
positive if any trajectory of the system starting at an initial 
state in the positive orthant R" remains forever in R". 

The following definitions of cooperative functions, homoge¬ 
neous functions and nondecreasing functions are widely used 
in the related works of investigating the stability of positive 
nonlinear systems with (or without time delays), like EU. 
Here we rewrite it to keep the self-integrity of this note. 

Definition 2: A continuous vector field / : R" —> R" 
which is continuously differentiable on R"\{0} is said to be 
cooperative if the Jacobian matrix df /dx is Metzler for all 
X e R+\{0}. 

Proposition 1: (see 1361) Let / ; R" —>■ R” be cooperative. 
For any two vectors x and y in R" with Xi = yi and x > y, 
we have ffix) > My). 

This Proposition can be obtained by the fact that for a 
cooperative function /, each (x) is a nondecreasing function 
of each Xj for j 7 ^ i. 

Definition 3: For an n-tuple r = (ri, • ■ • , ?-„) of positive 
real numbers and A > 0, the dilation map 5\ (x) is defined by 

When r = (I,!,-- - ,1), then dilation map is the standard 
dilation map. 


Definition 4: A vector field f{x),x G R", is homogeneous 
of degree p > 0 w.r.t. the dilation map S'^{x) if 

f{6lix)) = XW^ifix)) ( 1 ) 

Definition 5: A vector field g : R" —>■ R" is said to be 
nondecreasing on R", if for any x,y € R", and x > y, then 
9{x) > g{y). 


III. A SIMPLE BUT USEFUL TRANSFORM 


In this section, we first define a simple transform as: 


Zi 



( 2 ) 


where Vifi = I,-- - ,n is defined by the dilation map in 
Definition [ 3 ] For any function /, we can define a new function 
f{z), where 2 : = (zi, • • • , Zn), as: 


fiiz) 


MAP-XnP) _ MeI 

Ti-l — r.-l 7 

2 . 2 . 

0 , 


Zij^Q 
Zi = 0 


( 3 ) 


Next, we will present three good properties for this new 
function f{z) and prove them in the form of lemmas. 

Lemma 1: If f{x),x G R” is homogeneous of degree 
p > 0 w.r.t. the dilation map S^{x), then f{z),z G R" is 
homogeneous of degree p w.r.t. the standard dilation map. 

Proof: From the equation of ([T]i, one can get that: 




ffiX^^Xi,--- ,X‘^”'Xn) 


X^X^^Mx) 

An-i^p-i 


>^^SlMx)) 

X’^^-^zp-^ 


( 4 ) 


where /(O) = 0. The proof is completed. ■ 

Lemma 2: For the new function f{z) in Q, suppose the 
function / is cooperative, then for any two vectors z and w 
in R" with Zi = Wi and z > w, we have fi(z) > fpw). 

Proof: At first, from the definition of Q, if at least one of 
Zi,Wi is zero, the above claim is proved; otherwise, we have 


h{z) 



/*H 


fi{y) 


( 5 ) 


where Zi = x^’’' and Wi = . 

Therefore, if z > w and Zi = Wi, then x > y and Xi = yi, 
according to Proposition 1, one can get that fpx) > fpy) 
and zp~^ = 'wp~^. Combined with the equation (|5]), we have 
7i{z)>jpw). _ m 

Lemma 3: For the new function f{z) in Q, suppose the 
function / is nondecreasing and fi{x) = fl{xi), that is to 
say, there exists a function d{xi, • • • , Xi-i,Xi+i, • • • , Xn) > 
0, such that fi{x) > dxi. Then for any two vectors z and w 
in R” with z > w, we have ffiz) > ffw). 

Proof: In order to prove the conclusion, one can only 
need to check that dfPz)/dxj is nonnegative. 
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At first, from the condition that / is nondecreasing, one can 
get that 


In order to describe the convergence of the zero solution, 
we first define a norm as: 


dfijx) 

dxj 


> 0 , = 


,n. 


From Definition Q, we have fi{z) = fi{x)x[^ there 

fore, if i ^ j 

d7^{z) ^ a. 

Xz ^ yJ, 


dxj 


dxi 


if i = j S {1, • • • , n}, we have 


dxj dxj * 


=x. 


l/ri-2 




dxi 


Ti 


dfijx) ^ . l-n 


dxj ^ ~ r. 


+ ^/i(a;) > 


Therefore, one just need to prove that 
0. If not, simple calculations show that 

fi{x) < d{xi,- ■ ■ ,x^-l,Xi+l,■ • • ,Xn) ■ < dxi, 


which contradicts to the conditions fi{x) > dxi. 

Based on all above discussions, one can get that the function 
/j is nondecreasing. The proof is completed. ■ 


IV. Main Results 

At first, we give the network model for continuous-time 
positive nonlinear systems with unbounded time-varying delay 
as: 

/ iit) = f{x{t))+g{x{t-T{t)), t>0 

\ x{t) = t e (- 00 , 0 ] 

where Xi,i = I,-- - ,n is the state variable, g}(t) € 

C((—oo,0],M") is the vector-valued function specifying the 
initial condition of the system. r(f) is the time-varying delay 
with r(f) <t but unbounded. 

Assumption 1: Assume the following conditions hold: 

1 ) / is cooperative and g is nondecreasing on R", 

2 ) / and g are homogeneous of degree p w.r.t. the dilation 
map S'^{x). 

It isclear that these conditions are the generalization of 
linear positive systems, and /(O) = p(0) = 0, so 0 is a solution 
of system (|6]l. 

Using the above transform, let = Xi{t), then all 

Ti (t) 

+ giiz7 {t - r(f)), • • • , 4" (f - T(f)))] 

= ^-{7iiz{t))+g,iz{t-T{t)))] (7) 

where 7 i{.z{t)) = /*(x(f))/z[‘"\f) and - r(f))) = 

gi(x{t — T{t)))/zl' ^(f —r(f)). From Lemma [T] to Lemma 
12 we have the following property of /, and 'gp. 

1 ) / is cooperative and g is nondecreasing on M", 

2 ) / and g are homogeneous of degree p w.r.t. the standard 
dilation map. 


||2(0ll{L+oo} = . max ^Zi{t) (8) 

2 = 1 ,••• ,n 

Obviously, if ||2:(f)||{{,+oo} —> 0, then z{t) 0 and x{t) —>■ 0, 
as f ^ 0. 

Remark 1: In IZTl . a similar definition of weighted Zoo norm 
is defined by a vector u > 0 as 


\Xi 

max - 

l< 2 <n Vi 


It is clear that this norm has been introduced and used in 
discussing stability of neural networks successfully in ||29l, or 
even much earlier work ll30l and others. 

Definition 6: Suppose that g,{t) : IR.+ ^ IR+ is nondecreas¬ 
ing function satisfying /i(f) —>■ -boo as f —>• -boo. The positive 
nonlinear system (0 is said to achieve the global p-stability, 
if there exist constant scalars M > 0 and T > 0 (which may 
be related to initial values), such that for all f > T 

M 

ll^(i)llu.-roo} < ^, f = l, •••,«. (9) 

Remark 2: The concept of /r-stability was first proposed 
in ll2^ and for power-rate convergence in ll28l to discuss 
stability of dynamical systems with unbounded delays. It 
pointed out that between asymptotical stability and exponential 
stability, there exist various convergence concepts. It also 
shows that in case of bounded delay, the convergence rate is 
exponentially. And as the delay is increasing, the convergence 
rate is decreasing from exponential to asymptotical. 

Theorem 1: Suppose Assumption [T] holds. If there ex¬ 
ists a positive and nondecreasing function p,{t) satisfying 
lim p.{t) = -boo, and for j = 1, • ■ • , n 

t—>- + oo 


1 

r-i 


fM) 




— r(f)) ) 


+ lim 


flit) 


t->oo p,(ty-p 


< 0 . 


( 10 ) 


Then the zero solution of positive nonlinear system Q is 
globally p-stable, i.e., Zj(t) = 0{p~^{t)). 

Proof: From the inequality (flOl) . we can find a constant 
T such that for all t>T, 


IfYO , 

( tYo) ^ 



yp,{to - T{to)) J 

J 


-b 




< 0 . 


Moreover, for any x > 1, the following inequality holds: 


\fY0 

f tYo) \ 


(, * 

\p.{to - T{to))) 

6 J 


KtoY-p 


< 0 . 
( 11 ) 


Define the Lyapunov function as: 


V(f) = max{l,supU(f)}, (12) 

S<t 

where V{t) = p.{f) • ||-2:(f)||{{,+oo}- Then V(t) is a nondecreas¬ 
ing function, and V(t) < V{t). Now, we claim that V{t) is 
bounded. More precisely, we will prove that for all t > T, 

m = v{T). 
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In fact, for any time to > T, in case F(fo) < V(fo). then 
there exists a e-neighborhood of to, such that V{t) < V(fo) 
for t £ {to, to + e), i.e., V{t) is non-increasing at to- 

Instead, if V{to) = V(fo)j there exists an index j such that 

ti-{to)^J^Zj{to) = Vitoy,tiito)C^Zi{to) < V(fo),* ^ j, 

which mean that 

From the Assumption [T] using the cooperative and homo¬ 
geneous property of function /, one can get that 

= (.3) 

Moreover, since V{to — r(fo)) < V(fo) = V{to), therefore, 
- T{to))C^Zj{to - T{to)) < V{to), i.e.. 


Zj{to - r(fo)) < 


V{to) 




fl{to -T{to)) 

From the Assumption [T] using the nondecreasing and homo¬ 
geneous property of function p, one can get that 


gJz{to - T{to))) 


<9, 


V{to) 


^ \Kto - Tito)) 


< 




Now, differentiating I^(fo) by the upper-right Dini- 
derivative, and using the obtained results (fTTl i. (fT3l l and (fT4l i. 
we have 

,• Vito + A)-Vito), 

D+Vit)\t^to = lim sup---|t=to 

A->'0+ A 

=D+[nit) max 

i—l,--- ,n 

=D+[fiit)^-^z,it)]\t^to 
=Kto)^i^Zjito) + nito)^~^Zjito) 

7 ,:(tn'\ -I- ll.(f.n^ - 


=/i(fo)Cj ^Zjito) + nito) — [fjizito)) + gjizito - r(fo)))] 


gito) ^ 


(Vito)] 

\9ito)) 


p-i-i 


/,(e) 


-f 


Vito) 


gito-rito))) 


5,(0 


Vito) 

'vito)P 


( Pito) ] 


Pito)P 



\g.ito - Tito))) 

\ 


-f 


Ktp) 


< 0 , 


therefore, V(fo) is also non-increasing at to- 

In summary, V(<) = V(T) for all t >T, which implies that 
p(f)|| 2 ;(f)||{j^+oo} < V(T), i.e., the global /r-stability can be 
realized. The proof is completed. ■ 

If we define the norm for z(f) is as follows; 


ll^:(0ll{«.+oo,r*} = .max ^Ziit))^\ 


(15) 


then the corresponding results can be stated as: 


Corollary 1: Suppose Assumption [T] holds, if there ex¬ 
ists a positive and nondecreasing function /i(f) satisfying 
lim /i(f) = -bcx), and for j — - ■ ,n 

i—)-+oo 



( lim —— 

y t—>oo g,{^t 


9it) ] 

-Tit))) 



+ lim 

t—^OO 



< 0 . 


(16) 


Then the zero solution of positive nonlinear system © is 
globally p-stable, i.e., Zjit) = (f)). 

The proof is similar to that in Theorem [T] Here we omit 
the proof. 

Remark 3: By adjusting the parameter r*, the condition 
(fTbt can be satisfied. In lIZTl . r* = maxiry, while in ll23l . 
r* > max,; . 

In case p = 0 , the inequality (fTbl l can be rewritten as: 



( lim —— 

y t->oo 


9it) ] 
-rit))) 



lim M < 0 

t—)-oo 

(17) 


Except for the exponential stability and power-rate stability, 
we give the Log-stability and Log-Log stability (firstly given 
in 1291) criteria under this case. 

Corollary 2: (Log-stability) Suppose Assumption [T] holds 
with p = 0 , time delay T(f) < t — t/ \a.t, if 


/,(e)+5,(0 < 0 , (18) 

Then the zero solution of positive nonlinear system 0 is 
globally Log-stable, i.e., Zjit) = 0(ln(f -|- 1)“^). 

Proof: Choose the function p(f) = ln(< -f 1), and simple 
calculations show that 


p(f) 


lim 

t->oo gif! Inf) 


= 1 , and 


lim 

t^oo g,{t) 


Let r* = 1 , so the above inequality (fTTl i is equivalent to the 
inequality of (fTSl) . The proof is completed. ■ 

Corollary 3: (Log-Log stability) Suppose Assumption [T] 
holds with p = 0 time delay r(f) <t — t°‘,0<a< 1 , if in¬ 
equality (fTsT l holds, then the zero solution of positive nonlinear 
system 0 is Log —Log stable, i.e., Zjit) = 0(lnln(f-|-3) ^). 

Proof: Choose the function p(f) = lnln(f-|-3), and simple 
calculations show that 


limiiL 

t->oo /i(f“) 


= 1 , 


and lim = 0 . 
t -400 /i(f) 


Let r* = 1 , so the above inequality (fTTl i is equivalent to the 
inequality of (fTSl) . The proof is completed. ■ 


V. Numerical examples 
We choose the following example given by ET): 

iit) = fixit)) + gixit-rif)) (19) 


where 

—bx\ -f 2 a;ia ;2 
x1x2 — 4x2 


fixi,X2) = 
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and 


g{xi,X2) = 


X 1 X 2 

2xf 


In this system, it is easy to check that / is cooperative and 
g is nondecreasing; moreover, / and g are homogeneous of 
degree p = 2 w.r.t the dilation map 6^{x) with r = (1, 2). 

The difference of this example with that in lIZTll is that the 
time delay obeys the following condition; 

r(f) = t — t/ hit 

Let zi{t) = xi{t),Z 2 {t) = X 2 ^^(t). Then, after this 
transformation, the system ( fT^ turns to be 

zi{t) = —bzf + 2zi(f)z|(f) + zi{t — T{t))z 2 {t — r(f)) 

- r{t))] 

therefore. 



and 


f{Zl,Z2) 


—5zi + 2zi(f)z|(f) 

(z?(f)^|(f)-4z|(f))/z2(f) 


9{zi,Z2) 


zi(f)z|(f) 

2 z\(t)lz 2 {t) 


Obviously, these functions are homogeneous with degree 2 
w.r.t. the standard dilation map. 

If we choose /r(f) = ln(f + 1), ^ = (1,1) and r* = 2, then 


r mW r ln(f + l) 

lim ■ , , r = iim T—-7^ 

i-^oo/i(t/Ini) i-)-oo ln(i/Ini + 1) 


lim . . 

t^oo ^(i) 


AW 

1-1 


lim - 

t —>^00 X i 


= 0, 


( 20 ) 

( 21 ) 


and for j = 1,2, the left inequality of (fThl l in Corollary 1 
is —4 and —1 respectively. Therefore, the zero solution of 
positive nonlinear system ( fT^ can achieve the Log-stability, 
i.e., Zj{t) = 0(ln(f + 1)“^),_) = 1, 2, or equivalently. 


xi(f) = 0(ln(f-I-1) ^), X 2 {t) = 0{\n{t + 1) ^). (22) 


Choose the initial values as {xi{t),X 2 {t)) = (1,4) for t € 
(—00,0], then Fig. [T]shows that the trajectories of the positive 
nonlinear systems converge to zero. Moreover, even though 
the initial value of xi is larger than X 2 , however, from our 
above analysis, the variable X 2 (t) converges more quickly than 
xi (t). Fig. [T] also shows this phenomenon. Moreover, from 
the equation (|2|i, one can get that a better description of Log- 
stability is to plot the trajectories of ln(a;i(i)),* = lj2 w.r.t. 
ln(ln(f + 1)), see Fig. 2, from which one can see that our 
theoretical result does be a good approximation of the real 
convergence rate. 


VI. Conclusions 

In this note, the global /i-stability of positive nonlinear 
systems with unbounded time-varying delays is revisited. We 
propose a novel transform, by which the positive nonlinear 
can be transformed to a system with homogeneous nonlinear 
functions w.r.t. the standard dilation map. Then, within this 
framework, we analyze various /i-stability. Finally, a simple 
example is given to corroborate the effectiveness of the trans¬ 
form proposed in this note. 


Fig. 1. Trajectories of = 1,2 for positive nonlinear systems QD 

under unbounded time-varying delay r(£) = £ — £/ Int. 



Fig. 2. Trajectories of ln{xi{t)), £ = 1,2 w.r.t. ln(ln(£ + 1)) 
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